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We investigate the process of invasion percolation between two sites (injection and extraction sites) 
separated by a distance r in two-dimensional lattices of size L. Our results for the non-trapping 
invasion percolation model indicate that the statistics of the mass of invaded clusters is significantly 
dependent on the local occupation probability (pressure) p e at the extraction site. For p e — 0, we 
show that the mass distribution of invaded clusters P(M) follows a power-law P(M) ~ M~ a for 
intermediate values of the mass M, with an exponent a = 1.39 ± 0.03. When the local pressure 
is set to p e = p c , where p c corresponds to the site percolation threshold of the lattice topology, 
the distribution P(M) still displays a scaling region, but with an exponent a — 1.02 ± 0.03. This 
last behavior is consistent with previous results for the cluster statistics in standard percolation. 
In spite of these discrepancies, the results of our simulations indicate that the fractal dimension 
of the invaded cluster does not depends significantly on the local pressure p e and it is consistent 
with the fractal dimension values reported for standard invasion percolation. Finally, we perform 
extensive numerical simulations to determine the effect of the lattice borders on the statistics of 
the invaded clusters and also to characterize the self-organized critical behavior of the invasion 
percolation process. 



I. INTRODUCTION 

Multiphase flow phenomena in porous media are rel- 
evant in many problems of scientific and industrial im- 
portance, including extraction of oil and gas from under- 
ground reservoirs and transport of contaminants in soils 
and aquifers 0,0- From a theoretical point of view, these 
types of processes are very complex and hence difficult 
to be described in great detail. It is under this frame- 
work that pore network models have been used to repre- 
sent the porous media and concepts of percolation theory 
have been successfully applied to model slow flow of fluids 
through disordered pore spaces 0, 0, IE S. Q 13, 0, E3, EH • 
For example, the method frequently used for oil explo- 
ration consists in the injection of water or a miscible gas 
(carbon dioxide or methane) in one or more wells in the 
field, in order to displace and remove oil from the inter- 
stices of the porous rock. The extraction process persists 
until the interface of separation between the two fluids 
(water or gas and oil) reaches the extraction well where 
the injection fluid breaks through. At this moment, a de- 
crease in the oil production occurs. Because of econom- 
ical interests, it is important to determine, or at least 
estimate, the time when the oil production starts to de- 
cay. It is thus natural that companies of oil exploration 
demonstrate a great interest in the knowledge and de- 
velopment of this displacement process, specially in the 
investigation of injection dynamics, since the time asso- 
ciated with it is intrinsically related to the amount of oil 
that can be extracted from the underground. 

The production of oil in the petroleum field must be 
understood at the microscopic as well as at the macro- 
scopic level. During extraction, a decrease in pressure 



at the oil layer near the well usually leads to the move- 
ment of water and gas into this low-pressure zone. Ul- 
timately, when more distant oil cannot reach the well, 
water, gas, or both, are produced instead 0. In this 
case, the oil production is controlled at the macroscopic 
level, but eventually microscopic processes can also limit 
oil recovery. For example, capillary forces between water 
and rock can effectively close off small throats in a rock 
formation, so that oil cannot pass through them. The 
blockage at some pore throats can therefore prevent the 
extraction of oil in a large region of the oil field. 

It is important to determine the relevance of the in- 
volved forces in the displacement process. For this pur- 
pose, the capillary number is usually defined as 



where u is the flow velocity, /i is the viscosity of the 
displaced fluid and a is the interfacial tension between 
two fluids. This dimensionless parameter gives the re- 
lation between viscous and capillaries forces associated 
to the phenomenon. For high values of Ca, the viscous 
forces dominate the displacement dynamics. In the hy- 
pothesis of low values of Ca, capillary forces dominate 
and the dynamic of the displacement process is essen- 
tially determined at the level of pores, i.e., it is intrinsi- 
cally dependent on the local aspects of the geometry of 
the pore space. Under these conditions, heterogeneities 
of the porous media can be characterized by a random 
field, representing the spatial distribution of pore sizes, 
and the invasion percolation (IP) model can be applied 
to study the displacement of a fluid into a porous matrix. 

The IP model and its variants have been extensively 
used to simulate the process of displacement of a wet- 
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ting fluid through a porous medium by means of the 
injection of a nonwetting fluid with different viscosity. 
Such model has been very efficient when the injection 
process is quasi-static, i.e., in the regime of low veloci- 
ties. The phenomenon of displacement is then described 
by the growth of a cluster on a lattice, assuming that 
its border (perimeter) represents the interface of sepa- 
ration between the two fluids. Previous studies on this 
subject have shown that the cluster formed during the in- 
jection process is a fractal object with a fractal dimension 
close to that found in the traditional percolation model 
EHEII3 More recent studies have suggested that these 
dimensions are identical [l^. On account of its applica- 
bility, the IP model has been widely investigated in an 
attempt to determine its basic properties and to improve 
the predictions for problems of practical interest. 

In this work, the non-trapping invasion percolation 
(NTIP) model is applied in order to simulate the evo- 
lution of the interface between an invading and a defend- 
ing fluid in a porous medium between two sites (wells) 
separated by a distance r. We investigate the behavior 
of the mass distribution of the invaded cluster for dif- 
ferent values of the pressure in the extraction site. In 
this way, the amount of fluid extracted from the hetero- 
geneous medium is statistically quantified. The organi- 
zation of the paper is as follows. In Sec. 2 we briefly 
summarize the model used in the simulations. In Sec. 3, 
we present the results from computational simulations, 
while in Sec. 4 we conclude by discussing the significance 
of these results and possible applications. 

II. THE MODEL 

In the oil extraction process, when a wetting fluid (e.g., 
water) is injected slowly into a porous medium saturated 
with a nonwetting fluid (e.g., oil), capillary forces are the 
major driving forces. They determine the motion of the 
fluid when the injection process is carried out in the limit 
where Ca — > 0, i.e., the viscous effects are neglected in 
each pore, compared with the capillary effects. The IP 
model is useful in describing the basic features of this ex- 
traction process. In this situation, the pressure difference 
Ap between the two fluids across the meniscus follows the 
equation 

Ap=(^Jcos9, (2) 

where r p is the pore radius, 7 is the interfacial tension and 
9 is the angle between the interface of the two fluids and 
the pore wall. If we consider that 9 and 7 are constant, 
the advance of the separation interface will take place 
through regions of lower values of the capillary pressure, 
namely, those associated with large values of the local 
pore radius r p . 

In order to model the displacement of the interface 
between two fluids, we use the standard NTIP model in 
a two-dimensional square network. In this model, the 



displaced fluid is considered to be infinitely compressible 
and the injected fluid can penetrate through any region of 
the interface of separation into the fluid to be displaced. 
A regular lattice is adopted in the present study as an 
idealization of the porous medium and we assume that its 
sites represent elementary pore units. Following Eq. J5J, 
the microscopic features of the disordered porous medium 
are represented in this model through a random variable 
p ?» l/r p that gives the local accessibility of the pore 
space. 

Different from previous studies where the NTIP is ap- 
plied between two bars or from the center of the network 
to any of its edges, here the propagation front evolves 
from a site W\ (injection site) through the network un- 
til it reaches another site W2 (extraction site). These 
two sites W\ and W2 are located symmetrically between 
the limits of the network and separated by a distance r. 
Moreover, in our simulations the pressure at the extrac- 
tion site W2 is arbitrarily specified. We will show that 
the distribution of mass of the invaded clusters is signif- 
icantly dependent on this local value. For completeness, 
our NTIP algorithm is described as follows. Initially, a 
random value pi obtained from a uniform distribution 
limited in the range [0, 1] is distributed to each site of 
the network. In the initial step, the fluid to be displaced 
occupies all pores of the network, while the invading fluid 
is pushed through a single site, the site W\. We search 
among the neighboring sites of W\ the one which car- 
ries the smallest random number p. This site is then 
invaded, becoming part of the region occupied by the in- 
vading fluid, and the list of sites that are eligible to be 
invaded next is updated. The invasion process contin- 
ues until the site W2 is occupied, when the mass M of 
the invaded cluster is measured. For a fixed value of the 
local pressure at W2, p e , this process is repeated for sev- 
eral different network realizations and the distribution of 
mass of the invaded cluster P(M) is calculated. 



III. RESULTS AND DISCUSSION 

In Fig. 01 we illustrate some stages of the invasion pro- 
cess taking place in a lattice of size L = 256 and distance 
between sites r = 16 . The shades of gray indicate in 
which interval of time of the process a given site has been 
invaded. As can be seen from this typical realization of 
the model, the cluster resulting from the invasion process 
represents only a small fraction of the network. 

In what follows, we show how sensitive is the statistics 
of the mass of the invaded clusters to the pressure value 
p e imposed at the extraction site W2- Since the invasion 
finishes as soon as the extraction well is reached, we can 
interpret that the local pressure p e corresponds to a spe- 
cial criteria to stop the NTIP growth process, without 
any interference in the microscopic rules of the model. 
This suggests that the invaded clusters obtained for dis- 
tinct values of p e should still be in the same universality 
class observed for invasion percolation between two bor- 
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FIG. 1: Time evolution of the invasion process between the 
injection (W\) and extraction (W2) sites. The steps of inva- 
sion are identified in accordance with the number of iterations 
(time) executed during the simulation process. The colors in- 
dicate in which stage of the process the propagation front 
reached a given site. The light gray sites correspond to the 
range 1 < t < 1000, the dark gray sites to 1000 < t 5000, 
and the black ones correspond to i > 5000. The pressure 
value at the extraction well is p e — 0. 
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FIG. 3: Graph in log-log scale of the mass distribution of 
invaded clusters P(M) for the case in which the pressure at 
the extraction site is p e — 0, and a network size L — 512. 
Each curve correspond to a different value for the distance 
between the wells, namely r = 4 (plus), 8 (X) and 16 (stars). 
These distributions have been rescaled by its corresponding 
value Mo at the position of the bump. The straight line is 
the least-square fit to the data in the scaling region with the 
slope —1.39 ± 0.03. For better visualization the straight line 
has been slightly shifted upwards, i.e., the ordinate values 
have been multiplied by an arbitrary constant factor. The 
corresponding unsealed distributions are shown in inset. 
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FIG. 2: Log-log plot of the radius of gyration for L = 8192 
and r = 4. The circles correspond to p e = and the squares 
to p e = p c . These results have been produced with 10000 
realizations of NTIP clusters that do not touch the border of 
the lattice. The straight line (shifted downwards for better 
visualization) is the least-square fit to the data and the slope 
gives the fractal dimension df — 1.89 ± 0.01 . 



ders of a lattice. A direct evidence for this fact can be 
obtained through the calculation of the fractal dimension 
of the invaded cluster. We show in Fig.|2|the log-log plot 
of the mass of the invaded cluster against its radius of 
gyration for p e — and p e — p c . These results have been 
obtained for, L — 8192, r — 4 and 10000 realizations 
in which the growth process stopped before the invading 
cluster reached the border of the lattice. The best linear 
fit to both data sets gives approximately the same value 
for the fractal dimension df — 1.89 ± 0.01. This value 
is consistent with the value reported for the standard in- 
vasion percolation process, df ~ 1.89 [HISEIE] Our 
results therefore support the assumption that NTIP and 
SP are in the same universality class 

In Fig. |31 we show the log-log plot of the distribution 
P(M) of invaded masses for the casep e = 0, and different 
values of the "well" distance r. Clearly, all distributions 
displays power-law behavior P(M) ~ M~ a for interme- 
diate mass values. In addition, a lower cutoff of order r df 
is always present, corresponding to the minimum cluster 
size. The upper cutoff of order L 2 , on its turn, refers 
to the situation in which almost all sites of the lattice 
are invaded before the invasion front reaches the extrac- 
tion site. The effect of r on the distribution is to simply 
modify the range of the scaling region. The solid line 
in Fig. |31 corresponds to the best linear fit for the data 
corresponding to r = 4 in the scaling region, with the 
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FIG. 4: Log-log plot of P{M) for the case L = 256 and dif- 
ferent values for the distance between the wells, r — 4 (plus), 
r — 8 (X), and r = 16 (star). The pressure at the extraction 
site is equal to p e = p c — 0.5927. The solid line indicates the 
linear regression (for r = 4) with the exponent a — 1.00±0.03. 
For better visualization, it has been slightly shifted upwards. 



FIG. 6: Log-log plot of the distribution P(M) for p e = and 
L — 512. The squares include all clusters, the stars only the 
clusters that do not touch any of the edges, and the full circles 
only those that touch the edges of the lattice substrate. All 
distributions have been normalized by their integrals. 
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FIG. 5: Log-log plot of P(M) for the case r — 4 and different 
lattice sizes, L = 256 (plus), 512 (X), and 1024 (star). The 
pressure at the extraction site is equal to p e — p c ~ 0.5927. 
The solid line indicates the linear regression (for L = 1024) 
with the exponent a = 1.02 ± 0.03. 



exponent a = 1.39 ± 0.03. These curves also present a 
"bump" in the region of large values of mass. For com- 
parison, when we rescale the distributions by their cor- 
responding value at the maximum of this bump Mo, a 
collapse in the scaling region can also be observed. 

In Fig. 0| we show the distributions of mass P{M) for 
the case in which the pressure at the site W2 is arbi- 
trarily set to p e — p c ~ 0.5927, where p c corresponds to 



the site percolation threshold of the square network [jj. 
Again, the smaller the distance between the sites, larger 
is the range over which the power-law behavior holds. In 
addition, the distributions present the same lower cutoff 
of order r d f and the upper cutoff of order L 2 , where a 
bump can also be observed in the region of large values 
of mass ^l^j. We show in Fig. [3] the effect of the sys- 
tem size on the distributions of mass for the case r = 4 
and p e — Pc- The results shown in Figs. 0] and [S] are also 
consistent with the existence of a power-law regime for 
intermediate values of mass, P(M) ~ M~ a , but with an 
exponent a — 1.02 ± 0.03. The less negative exponent 
in the power-law distribution for the cases p e — p c when 
compared to p e < p c means that, in the later case, larger 
clusters appear more frequently. 

As already mentioned, just before the extraction site 
is invaded, all sites belonging to the invasion front have 
their pressure values larger than p e . Therefore, in the 
case where p e = p c , the invaded cluster is similar to 
a cluster generated using the standard percolation (SP) 
model at the critical point. In the SP model the distri- 
bution of cluster sizes follows a power-law, P(S) ~ S~ T , 
where the exponent r m 2.055 is known as Fisher expo- 
nent 0. Assuming that the IP and SP models are com- 
patible, note that, in our model the invaded cluster would 
be the same if any of its sites were chosen as the injec- 
tion point. As a consequence, the probability of finding 
an invaded cluster of size M is the product of two factors, 
namely, the probability with which such cluster appears 
in a percolation lattice at the critical point, and the size 
of the cluster. Thus, we should expect a — t— 1 sa 1.055, 
a value that is in good agreement with our numerical re- 
sults. 
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FIG. 7: Log-log plot of the distribution P(M) of the invaded 
cluster p e — p c and L — 512. The squares include all clus- 
ters, the stars only the clusters that do not touch any of the 
edges, and the full circles only those that touch the edges of 
the network. All distributions have been normalized by their 
integrals. 



FIG. 8: The fraction \ of invading clusters that do not touch 
the edges of the network as a function of the pressure value 
at the extraction site for system size L — 256 and distance 
between sites r = 8. 



The bumps present in the distributions shown in 
FigS-EJand^are due to aggregates that grow long enough 
to reach the border of the lattice before reaching the ex- 
traction well. This is made clear when we distinguish 
these clusters in our simulation data from those that stop 
before reaching the borders. The distributions shown in 
Fig. Eland Fig. reveal that the bumps arc mainly origi- 
nated by the cluster masses subjected to size limitations 
(full circles), i.e., those clusters that touch the lattice bor- 
ders. Note that, even after touching the border, the clus- 
ter continues to grow until the extraction well is invaded. 
Thus, those clusters that contribute to the bump in the 
region of large M have a fractal dimension dt = 2 and 
therefore do not possess the critical properties observed 
in percolation clusters. The statistics of these larger clus- 
ters can be investigated if we compute the fraction \ of 
invaded clusters that do not touch the edges of the sim- 
ulation box for different values of the pressure p e . As 
shown in Fig. [8] for a given value of L and r, \ remains 
approximately constant up to a value of p « p c where 
it experiences a sharp decrease and reaches x — 0- F° r 
p e < Pc, the presence of these large clusters is an artifact 
of the finite size of the lattice and becomes less relevant as 
the lattice size grows. In the results which follow we will 
only consider the clusters that do not touch the borders 
when computing the distributions of mass. 

To better understand how the statistics of the mass dis- 
tribution changes with the value of the pressure at the ex- 
traction site, we perform simulations with different values 
of p e = 0.0, 0.2, 0.4, 0.55, p c (~ 0.5927), and 0.6. From 
the results presented in Fig. EH we observe that the dis- 
tribution P{M) remains practically unchanged (see the 



inset) with a scaling regime at intermediate values of M 
that persists at least up to p e — 0.4. For p e values larger 
than this, the distribution P(M) changes gradually (see 
the curve for p e = 0.55) until another power-law is effec- 
tively established at p e — p c . For p e > p c , most of the 
generated clusters are subjected to border effects and the 
scaling behavior is again suppressed. 

The behavior of the distributions as p e approaches p c 
is shown in Fig. ^| For large values of M, all curves 
converge to the same behavior, which coincides with the 
power-law form obtained for the case p e = 0. For small 
values of M, the curves display a power-law decay with an 
exponent that depends on the value of p ei but approaches 
the value a = 1.055, expected in the limit p e = p c . 
One can understand this behavior by making a parallel 
with the way in which the standard percolation model 
approaches the critical point. As p — > p c in SP, the char- 
acteristic length of the clusters £ diverges as (p c — p) -ly , 
with v = 4/3. Thus, we should observe the critical be- 
havior for values of M bellow the characteristic mass of a 
cluster obtained with p — p e , and the off-critical behavior 
for values of M above this characteristic mass. It means 
that the onset of the crossover between the two behav- 
iors should scale as {p c — Pe) -7 , with 7 = vdf. Assuming 
that our invasion process is consistent with SP, we ex- 
pect 7 w 2.52. This ansatz is supported by the scaled 
curves shown in the inset of Fig. ^| As can be seen, for 
Pe <C Pc, the crossover moves to values of mass smaller 
than the cutoff region in r f , and only the off-critical 
regime is present. On the other hand, for p e > p c , most 
of the generated clusters reach the border of the lattice, 
leaving the critical regime, and the scaling behavior is 
again suppressed. 

It is possible to understand why we obtain distinct 
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distributions in the critical (p e — p c ) and off-critical 
(Pe < Pc) regimes by noting that our invasion percola- 
tion process can be divided in the following two sequen- 
tial phases. In the first phase, the invaded cluster grows 
from the injection site and reaches an immediate neigh- 
borhood of the extraction site. The second phase includes 
all additional invading steps up to the point in which the 
extraction site itself is taken, and the invasion process is 
then terminated. Thus, we can say that the total invaded 
mass is the sum of two variables, M = Mi + M 2 , where 
Mi and M 2 are the numbers of invasion steps to accom- 
plish the first and second phases, respectively. Since the 
masses Mi and Mi have distinct distributions, namely, 
Pi (Mi) and P%(M?), respectively, the distribution with 
the longer tail between them, should determine the be- 
havior of the tail of P{M). Note that the distribution 
Pi (Mi) does not depend on the value of p e . On other 
hand, the number of additional steps to invade the ex- 
traction well clearly depends on p e and so does P2(M 2 ). 
In the particular case of p e = 0, it follows that M 2 = 1. 
We can then write that P(M) = P X (M - 1) ~ M~ a , 
which shows that the exponent a — 1.39 comes from 
the distribution of steps necessary to complete phase 1. 
When p e approaches p c , the contribution of M 2 to the 
distribution P(M) starts to become relevant. In the crit- 
ical condition p e — p c , the tail of this distribution is 
dominated by the number M 2 of steps to invade the ex- 
traction well. Although we lack information about the 
complete form of P 2 (M 2 ), we know that the necessary 
condition to invade the extraction well is that all sites 
belonging to the perimeter of the invaded cluster have 
pressures p > p e . Following percolation theory when 
Pe — Pc, this condition leads to a power-law distribution, 
P(M) ~ M~ a , where a = 1.055. As a consequence, 
we can deduct that, to be dominant, the tail of P 2 (Af 2 ) 
should behave in the same way. 

The invasion percolation is broadly recognized as the 
most convincing model for a natural phenomenon that 
displays self-organized critical (SOC) behavior [l^. In- 
deed, this characteristic is manifested, for example, in the 
distribution of the random numbers P(p) corresponding 
to those sites that have been occupied during the inva- 
sion process. Accordingly, the fact that P(p) displays 
a sharp cut-off on the right side at the critical percola- 
tion point p c represents a sound signature of the SOC 
behavior, which dictates the dynamics of the displace- 
ment process. In Fig. ^] we plot the distribution P(p) 
computed from 1000 realizations of size L = 512 and dis- 
tance r = 4, and for values of p e = 0, 0.61, 0.63 and 0.65 
as well as for p e = p c . For p e — and p e = p c , after the 
completion of the invasion process, both distributions dis- 
play a clear cutoff at p ~ p c , in agreement with previous 
simulations of the standard invasion percolation model 
pbl l2lj. The inset in Fig. 1111 shows that the transition 
is clearly sharper for p e = p c than for p e — 0. In the 
cases in which p e > p c , we observe that the cutoff of the 
distribution P(p) shifts rightwards to p — p e . Since we 
are only interested in those clusters that do not touch 
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FIG. 9: Log-log plot of the mass distribution of the invaded 
cluster for different values of the pressure at the extraction 
site. p e = 0.4 (circles), p e = 0.55 (full circles), p e = p c 
(triangles) and p e — 0.6 (full triangles). The inset shows the 
mass distribution of invaded cluster for pressure values p e = 
0.0 (circles), p e — 0.2 (triangles) and p e = 0.4 (diamond). 
The system size is L — 256 and the distance between the 
injection and extraction sites is r — 8. 



the borders, for p e > p Cl a very large number of trials 
of the invasion process is necessary to generate the num- 
ber of realizations required to produce good statistics (see 
Fig.|SJ. Once more, if we argue that p e can be interpreted 
as a stopping criteria for the invasion process, the non- 
invaded sites that are neighbors of the perimeter sites of 
clusters generated under such constraint should all have 
local probabilities that obey the inequality p > p e . As a 
consequence, for p e > p c , the invaded clusters are limited 
by a "high energy" border of non-invaded sites. 



IV. CONCLUSIONS 

In summary, we investigate here the statistics of the 
mass of the invasion percolation clusters generated be- 
tween two sites separated by a distance r in a two- 
dimensional network, and for different values of the local 
probability at the extraction site, p e . The results of our 
simulations indicate that the distribution of mass of the 
invaded cluster P(M) is highly sensitive to p e . More 
important, P(M) displays power-law behavior in the fol- 
lowing two different conditions: (i) an off-critical regime 
for sufficiently low values of p e , and {ii) a critical regime 
for p e = p c - Interestingly, even though the scaling expo- 
nents calculated in both cases are substantially different, 
namely a — 1.39 ± 0.03 and 1.02 ± 0.03 for cases (i) and 
(ii), respectively, our results for the fractal dimension of 
the invaded clusters show that they should belong to the 
same universality class. For p e approaching p c , we found 
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that the distribution displays a crossover from the crit- 
ical to the off-critical regime and that the onset of the 
crossover scales as £ df , where £ is the correlation length, 
£ ~ (Pc — Ve) v - Finally, we have shown that the SOC 
behavior of the invasion percolation process is rather ro- 
bust, even in the situation in which p e > p c . In this 
case, although rarely found, those clusters whose growth 
has not been affected by the borders of the NTIP system 
(named here "critical clusters'' ) present a very peculiar 
property, namely, that all the non-invaded sites at their 
perimeters have an anomalously high local probability. 



FIG. 10: Log-log plot of the mass distribution of the invaded 
cluster for values of the pressure at the extraction site ap- 
proaching the critical point p e — p c for systems size L = 8192. 
The thin lines are included as guide to eyes and have the 
slopes expected in the following limiting cases: (i) —1.05 for 
the critical case p e = p c and (ii) —1.39 observed in p e = 0. In 
the inset, the curves are rescaled by the characteristic mass 
£ d f , where £ is the correlation coefficient which scales with 
the pressure at the extraction well as £ ~ (p c — p E )~ v ■ The 
exponents d/ and v are assumed to be the same as in the the 
standard percolation model, namely, df = 1.89 and v = 4/3. 
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FIG. 11: The probability distribution P(p) of invaded sites for 
different values of the pressure in the extraction site, p e = p c 
(circles), 0.61 (squares) and 0.63 (diamonds) and 0.65 (trian- 
gles) for L = 512 and r = 4. The inset shows the probability 
distribution P(p) forp e = 0.0 (dotted line) andp c (solid line). 



